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Abstract 

Twenty years of extensive research has yet to produce a general consensus on the 
origin of high temperature superconductivity (HTS). However, several generic char- 
acteristics of the cuprate superconductors have emerged as the essential ingredients 
of and/or constraints on any viable microscopic model of HTS. Besides a T c of order 
100 K, the most prominent on the list include a d-wave superconducting gap with Fermi 
liquid nodal excitations, a d-wave pseudogap with the characteristic temperature scale 
T*, an anomalous doping-dependent oxygen isotope shift, nanometer-scale gap inho- 
mogeneity, etc.. The key role of planar oxygen vibrations implied by the isotope shift 
and other evidence, in the context of CuC>2 plane symmetry and charge constraints 
from the strong intra-3d Coulomb repulsion U, enforces an anharmonic mechanism in 
which the oxygen vibrational amplitude modulates the strength of the in-plane Cu-Cu 
bond. We show, within a Fermi liquid framework, that this mechanism can lead to 
strong d-wave pairing and to a natural explanation of the salient features of HTS. 
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Great strides towards understanding HTS have been made during the last twenty years of 
intense experimental and theoretical research on cuprate superconductors 

HQ 

. However, 

there is no general consensus on the origin of HTS. Especially there is no report on a specific 
microscopic pairing mechanism that is capable of encompassing the complex phenomenology 
of the superconducting and normal states consistently in one theoretical model. 

Among the highly unconventional properties of the cuprate superconductors are the d- 
wave symmetry of the superconducting gap jjj and the presence of a pseudogap also with 
<i-wave symmetry Q] . Theoretical insight is provided by study of the low-energy excitations 
around the node in the ci-wave gap, where disparate low-temperature experiments including 
specific heat |5[ , transport 



,1, 3, Q, Q, 



and penetration depth |12| , together with 
angle-resolved photoemission spectroscopy (ARPES) 13] can all be interpreted in terms of 
a Fe™ li q oid deseripUon. This suggests the plausibihty of a broadly BCS f_k g 
of which an important consequence is that the large on-site Coulomb repulsion U jlj] does 
not enter into ci-wave pairing to first order. 

In conventional superconductivity the isotope exponent a has been key in signalling the 
role of phonons in the pairing mechanism. HTS exhibits a universal, anomalous doping- 



dependent isotope shift 



3, m 



which shows that phonons are again playing a key role, but 



in some unconventional manner, and any viable microscopic model must be able to account 
for this. A further challenge to theory is the recently observed spatial inhomogeneity of the 
gap at the nanometer scale 

Recently interest in the role of electron-phonon coupling has been rekindled by mounting 
experimental support, which, in addition to t 



: re isotope shift, includes electronic Raman 
, inelastic neutron scattering j^J, x-ray 
absorption fine structure |25|, low-temperature STM/STS and isotope effect in pen- 



scattering, 
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etration depth measurements |2 
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Extensive theoretical studies |2£ 



Q El 



have shed light on the relevance of electron-phonon interaction in 
understanding HTS. 

Here, in re-examining phonon coupling, we find a novel microscopic pairing mechanism 
which can indeed explain the essence of HTS phenomena: high T c 's, <i-wave pairing, d- 
symmetry pseudogap, anomalous isotope shift, and nanoscale gap inhomogeneity. Although 
the effects of U are manifested |[1J], here we shall focus on interactions directly responsible 
for pairing, allowing for U indirectly by decoupling from on-site charge fluctuations. We 
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FIG. 1: (a) The unit cell in the Cu02 plane, with Cu atoms (yellow) and O atoms (red), showing x, 
y, and z vibrational modes (arrows), (b) The Cu 3d x 2_ y 2 and O 2p x orbitals, illustrating effect of 
an O z-displacement (green arrows), positive sense (top panel), and negative sense (bottom panel). 

start by considering how oxygen vibrations can modify electronic motion in the Cu02 plane, 
the universal active component of HTS materials. The Cu02 plane is a square lattice of 
divalent Cu ions with oxygen ions located at the centers of the Cu-Cu bond (Fig. la). Only 
the highest-lying 3d x 2_ y 2 orbital (Fig. lb) plays an active role, crystal field effects demoting 
the other Cu 3d 9 orbitals to core-like status. 

Let us look at a single Cu-Cu bond between 3d x 2_ y 2 orbitals located on atoms labeled 1 
and 2 (Fig. 1). The kinetic energy associated with the unperturbed direct bond 1-2 comes, 
within a one-band model, from the nearest-neighbor hopping matrix element t 

h e u = ~t ( c t« c ^ + c %c c rA > (!) 

where the cf a (cj i(7 ) are Fermion creation (destruction) operators for Cu site i and spin a. 
The intersite hopping actually occurs via superexchange, i.e. wavefunction overlap between 
the 3d x 2_ y 2 orbital of Cul, the intrabond oxygen 2p x longitudinal with the bond, and the 
3d x 2_ y 2 orbital of Cu2 (Fig. lb), sensitizing it to the local oxygen vibrational degrees of 
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freedom (Fig. la). Consider the effect, sketched in Fig. lb, of, for example, the out-of- 
plane oxygen displacement z, which is to reduce the overlap between both 3d x 2_ y 2 orbitals 
and the oxygen 2p x orbital, thus reducing the effective coupling t. But due to the local 
centrosymmetry of the O-site, this t-reduction effect is the same irrespective of the sign of 
the displacement z, so that its lowest order expression is as z 2 . Hence the electron-vibrator 
term in the Hamiltonian must have the unusual second order form of coupling 

h ii = -/ E + c t c ^) > ( 2 ) 

cr 

where v is the coupling strength. This electron-vibrator coupling causes the Cu-Cu bond 
strength t to fluctuate with the oxygen vibrator square amplitude z 2 - hence the description 
Fluctuating Bond Model (FBM). The two other oxygen vibrational modes, x and y, can 
also couple to the bond strength in a similar manner. In addition to coupling to the bond, 
the vibrational mode x, longitudinal to the bond, can also displace charge onto the Cul 
and Cu2 sites; however since charge accumulation on the Cu sites is resisted by the large 
intrasite Coulomb interaction U, we do not believe that on-site charge displacement can be 
important, and ignore it in our model at the present maximally simplified stage. 

The oxygen vibrational degree of freedom needs to embody an anharmonic potential 
term in a theory with nonlinear electron- vibrator coupling, e.g. for the z-mode the vibrator 
Hamiltonian should have the form 

where a fourth-order potential with coefficient w (again, centrosymmetry is used to eliminate 
a cubic term) has been included, in addition to the conventional harmonic terms in which 
p z is oxygen momentum, m effective mass, and uq the harmonic vibrator frequency. The 
potential may be a flattened well or a double well, depending on the sign of Uq (see Fig. 
2a). A quartic anharmonic potential is not unique in perovskite physics, a similar potential 
plays a key role in the theory of ferroelectricity 37 1 . 

In the following the FBM Pairing Interaction will be derived systematically. First, we 
give a non-rigorous argument leading to a physical picture of the pairing by completing the 
square in the highest-order interaction. Introducing the compact notation as the bond 
order operator, describing the strength of the Cul-Cu2 electronic coupling, 

X 12 = E i C taC2,a + < CT C V ) , (4) 



FIG. 2: Left (a): Bare oxygen anharmonic potential (3) full (dashed) curves with positive (negative) 
harmonic coefficient uiloq. Right (b): Binding of two quasiparticles (e.g. holes h) in a Cu-Cu bond, 
as given by Eq. (5). 



and completing the square (classically valid in the Fig 2a double well case when —muj^ » 
v(X 12 )) 

-z* + V -^X 12 = ™ (V + ^X 12 ) 2 - *Xt 2 , (5) 
8 2 8 V w J 2 12 ' V ; 

we find a new, attractive electron-electron (or hole-hole) interaction in the intra-bond 
channel, of strength K, where K = v 2 /w. K is the key interaction in the FBM. Because 
the interaction K acts via the bond operator X\ 2 , it is seen from the form Xf 2 to lie in the 
intra-bond channel (see Fig. 2b) as opposed to being on-site, and with the |/-bond included 
the Fourier transform of the interaction leads to the d-wave factor cosq x — cosq y , resulting 
in pairing and pseudogap phenomena of d-wave symmetry. There is experimental evidence 
for linking the HTS pairing mechanism with the oxygen ^-vibrational mode: if the vibrator 
is localized in the z-direction by large static Cu0 2 plane buckling 3^, T c is reduced and 
can go to zero. Hence a key role in pairing must indeed be played by the oxygen degree of 
freedom. 

We are now ready to write down the complete FBM Hamiltonian, as a sum of electronic, 
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vibrator, and coupling terms: 

jjFBM = H e + H v + H ev (q) 

Here the electronic term includes hopping over longer ranges than the nearest neighbor 
hopping considered in EqtJTJ): 

H e = -\Y,t(\-])cf a c^ (7) 

where i denotes the 3d x 2_ y 2 orbital on lattice site i = (i x ,i y ) in the 2D square lattice of 
Cu ions. The strongest interaction is the nearest neighbor hopping integral £(±1,0) = 
t(0, ±1) = t, followed by the next-nearest neighbor interaction i(±l, ±1) = t', and then the 
3rd-nearest neighbor interaction £(±2, 0) = t(0, ±2) = t" . The band eigenvalues of (J7J) 
are = — 2t(cos k x + cos k y ) —At' cos k x cos ky — 2£"(cos2/l£ -\- cos 2/c^) in units where lattice 
const ant = 1. 

In the vibrational term 



1 2 mul 2 W , 2 y 

o Pi.a ~l ^ X i,a 7, \ X i,a) 



(8) 



the bonds are relabelled for greater notational convenience, in terms of a Cu site i and a 
direction a = x or y away from i in the positive axis direction. The vibrational modes are 
approximated as local (Einstein) and isotropic. The notation p 2 , x 2 implies Y^gPl, ^ s ij 
where s is polarization (s = transverse to plane, longitudinal to bond, or in plane transverse 
to bond). The mode degeneracy is n. 

The bond order operators are defined using the same bond notation 



X i>a = [ C ta C i+Z a ,° + CRx^a- , (9) 



where x Q is a unit vector along the direction a. In terms of these, the coupling Hamiltonian 
is 

).. Inn < ' ' 
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The prefactor includes a spin degeneracy n s . 

To solve the FBM we shall in this paper use the standard weak-coupling approach, based 
on an electron gas as the unperturbed system. To develop a perturbation expansion in the 
absence of the Migdal theorem in HTS Q], we adopt the 1/N expansion technique, where 




FIG. 3: (a) The interaction between two pairing quasiparticles k and — k, exchanging momentum 
q and frequency u> n , is a product of <i-wave form factors C^k+q (yellow circles) and potential 
propagator V (q, cj n )(b) First approximation to V based on v 2 x two-boson propagator. Dashed 
lines represent single vibrator (boson) propagator, red circles electron-boson interaction v. (c) The 
.eading-TV self-energy corrections to V (see analogous corrections in the heavy fermion problem 



41 1) come from boson-boson interaction w and v 2 x Response Function Rdd- Full lines represent 
fermion quasiparticle propagator. Green square represents boson-boson interaction w, fermion lens 
is response Function Rdd- 

N is degeneracy, e.g. orbital or spin degeneracy. The 1/N expansion works well e.g. for 
the Kondo problem 3^|, the results remaining physical down to spin degeneracy N = 2. 
Here we systematically co-expand in the inverse of the mode degeneracy n and the spin 
degeneracy n s using a path integral approach (see supplementary material), meaning 
by expressions such as "1/N" the joint orders 1/n and l/n s . 

The interaction which scatters a pair (k, — k) to (k + q, — k — q) (Fig. 3) differs from 
the usual single phonon propagator structure in BCS, because there is no single electron- 
phonon interaction term in the FBM Eq's ([BlllO)). The structure is instead a two-boson 
propagator (Fig. 3b) with self-energy insertions, the leading-iV self energy involving the 
boson-boson interaction w and v 2 x a fermion response function (Fig. 3c). Because the 
response function can go from bond a = x or y to bond (3 = x or y it is a 2 x 2 matrix. 
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FIG. 4: The Fluctuating Bond field after projecting out s-fluctuations has d-symmetry, with (left 
panel) O-amplitude (dumbells) large in y-bonds, small in x-bonds, the reverse in right panel. Panels 
also illustrate opposite phases of the static d-wave CDW (18). 

The matrix can be simplified by noticing the dramatic divergence produced by the self- 
energy insertions at small transferred wavevector q. This small-g divergence comes from 
fluctuations having o?-wave symmetry, the s-wave part being nonsingular and uninteresting. 
Therefore, our simplifying procedure is to reduce the 2x2 propagator matrix to a scalar 
by projecting out s-symmetry charge fluctuations and retaining only the (^-symmetry ones, 
the projection procedure being implemented in path integral formalism (see supplementary 
material). The vibrator amplitude/X-operator fluctuations around each Cu site now form a 
d x 2_ y 2-Y\ke pattern (see Fig. 4). This selection of the <i-channel has the additional physical 
merit that charge flow into and out of each Cu site is balanced (Fig. 4), i.e. there is zero 
net site charge accumulation, compatible with the accepted large Coulomb repulsion U on 
each site, which inhibits such charge fluctuations. 

The FBM pairing propagator r(k, q, n) (Fig. 3) for scattering a pair from (k, — k) to 
(k + q, — k — q) is now given (for a full derivation see supplementary material) by the Fig. 
3 graphs as a product of <i-wave form factors and a 2-phonon potential propagator V(q, n) 




(11) 
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where the <i-type form factor (coming from the form factor of the bond operator X) is 

Ck,k' = 2 [ cos + cos (K) - cos (K) ~ cos (K) ] > ( 12 ) 

and the potential propagator (dropping a zero-frequency term, only significant at tempera- 
tures above those of interest) is 

Vfq n) = -^KulU 

^ J ^ + 4^ + 4^ /q [I-^(q,n)]- 1 " J 

Here o; n = l^nk^T is the Matsubara frequency (Fourier component with respect to imag- 
inary time), T =temperature, ks =Boltzmann's constant, and we introduce definitions of 
the total vibrator frequency 57, the mean field harmonic vibrator frequency Uh, and the 
anharmonic component of the vibrator frequency u a 

tf = ul + ul, (14) 



u 2 + - 



m 

In the denominator of (jlHj) we can identify (in square brackets) the two terms in the self 
energy (Fig. 3c), the positive one coming from w, and the negative one coming from v 2 Rdd- 
The other two terms are the inverse of the two-boson propagator (Fig. 3b) with characteristic 
frequency 257 (a frequency up-shifted by the single-boson self energy). The form factor / q is 
defined by / q = ~ [cos 2 (q x /2) + cos 2 (q y /2)] , and the dd Response Function R dd is defined 
in the normal state as 

D ( \ /( e k) — /(fk+q) c c 

^ iu n + e k - e k+q 

where /(e k ) is the Fermi function. The dd Response Function (RF) is a generalization of 
the density-density RF. The presence of the ^-factors shows that the RF can be interpreted 
as the (i-symmetry density response to a ci-symmetry perturbation. At long wavelengths, 
the RF is decoupled from the on-site fluctuations of s-symmetry, hence it is decoupled from 
the strong on-site Coulomb repulsion U. 

The modified 2-phonon interaction V(q, n) can diverge at low frequency and small q, at 
significant values of the interaction K, due to the largeness of Rdd when the Fermi level lies 
at the energy e SP of the saddle points (SP) (or "antinodal points") at k = (n,0) and (0, it) 



in the band structure, signaling the emergence of a 2-phonon bound state. The same SP 
effect causes a peak, the van Hove singularity, in the density of states (DOS) at energy e$p- 
In this situation the superconducting gap acts, through controlling the magnitude of R^d, 
to regularize the divergence. 

Let us now look at the leading- N gap equation at T c 

A (k,n) = ~TJ2 £U' V ( k - k '> n ~ n ') G * ( k >0 A ( k >') > ( 16 ) 

k',n' 

where A (k,n) is the gap and G2 (k,n) = [y^ + e|) \ 

In standard BCS theory, G2 gives rise to a log divergence in T, which results in a solution 
to the gap equation with the well-known standard BCS formula for T c exponential in the 
inverse coupling constant. In contrast, in the FBM there is the stronger divergence in 
V(q, m), which tends to peg T c at the temperature where the divergence disappears. Since 
this temperature is defined by the denominator in V(q, m), in which no degeneracy factor n s 
appears, it is a leading- N formula. Hence we are justified in neglecting the second Eliashberg 
equation involving electronic mass renormalization (Z-factor) as a 1/N correction. 

Solving the gap equation, we always find a d x 2_ y 2 -wave gap A (k,n) due to the k , factors 
- originating in the FBM pairing interaction being localized on a bond (Fig. 2b) - and to 
the largeness of the V(q, n) at small q. In Fig. 5a we present some results for T c and oxygen 
isotope shift (the technique used was to solve at finite gap using Fast Fourier Transform 
(FFT) techniques and extrapolate to zero gap) as a function of doping. The results 
show the standard hump in T c as function of doping, and a very dramatic minimum in the 
isotope shift going down to almost zero around the T c -maximum, while going up to values 
above the BCS a = 0.5 on the underdoped side. In Fig. 5b we present experimental results 
for the "universal" isotope shift behavior found for several materials Q], along with the 
widely used "universal" empirical formula for T c 43]. It is seen that there is a remarkable 
semiquantitative agreement between theory and experiment, which is especially remarkable 
as the isotope shift in BCS theory is quite robust. 

The explanation for the hump in T c as a function of doping in the FBM, is that the 
density of states and the RF peak around the point where the Fermi energy coincides with 
the band structure energy esp at the SP. At this point a zero in the denominator of the 
pairing interaction V(q, m) is more readily obtained, while high DOS always favors pairing. 
The phonon parameters in Fig. 5 have been selected so that the phonon frequency is 
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FIG. 5: (a) Left panel. Transition temperature T c and oxygen isotope shift a vs. doping p, 
relative to doping po at max. T c , from gap equation (16). Parameters t = 0.25 eV, t' = —0.06 eV, 
t" = 0.0325 eV, K = 0.48 eV, uj a = 0.05 eV, u>h = 0.015 eV. Phonon frequencies in this range are 
reported by ref's ^^|-j3](b) Right panel. Experimental transition temperature and isotope shift. 

n 

Full curve, empirical T^-doping relation T c — T c max 

(1 - 82.6 (p - po) ) 43]. Points are oxygen 
isotope shift measurements for YBCO-based and Bi-2212 cuprate superconductors . 



mainly of anharmonic origin, i.e. a mainly quartic potential (see Fig. 2a), when the zero 
in the denominator of V(q, n) is mass-independent (u a factors out), hence the extremely 
low isotope shift. As the Fermi energy moves away from esp, the DOS drops and the FBM 
divergence in the pairing interaction tends to disappear, with a resumption of more normal 
BCS isotope shift. 

As well as driving pairing, the coupling K can also produce static distortions. The 2- 
phonon interaction w can be renormalized to w by summing a somewhat more extended set of 
leading- N Feynman diagrams than those in Fig. 2. w is found to contain the proportionality 
factor 

w ~ ~ - KR dd (q, n), (17) 

leading to a zero in the long- wavelength w at a temperature T m f. Below T m j there is 
symmetry-breaking in the system, which can be likened to the presence of an Ising pseodospin 
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in each unit cell. The symmetry breaking can be described in real space as a local splitting 
of the nearest-neighbor hopping integrals t, t x ^ t y , and in k-space as a splitting of the 
saddle point energies at the fc-points X and Y, ^ €g P . 

An ordered pseudospin structure has been sought in the form of a ID d-wave CDW j^J , 
with the Ansatz 

9 2Jnn s . .. , , 

u\ = ^ r ^x Q cos(Q.i), (18) 

where 

2 — ^ ( 2 I 2 2 2\ / ln \ 

M i — 2 V X i+x/2 + x i-x/2 ~ 2A+y/2 — Vi-y/2) > l iJ J 

and Q = (Q x , Qy) is the CDW wavevector. The c?-wave nature of the CDW is seen in that 
the expression (fTTJj) corresponds to a static distortion of the type represented in Fig. 4. 

The propagator G(k, n), in a standard approximation where the self-energy is second 
order in xq, is given by 

G(M) = *,-»-W(M) - (20) 



with 

i T A 2 , . „ A 2 , „ 1 

(21) 



$(k,n) = ^ 



^k,k+Q 'Ck.k-Q 



_Wn — Ck+Q ^^n — Ck-Q. 

By calculating the free energy, the wavevector Q was determined by energy minimization. 
The direction of Q is found to lie along the x- or y- axis, Q = (Q, 0), or = (0, Q). The CDW 
is found to be stable - relative to the uniformly-polarized (but presumably disordered) state 
- below a temperature Tcdw, which is plotted in Fig. 6, a stability mostly confined to the 
underdoped side. Q in this region closely tracks the nesting wavevector between the two 
pieces of Fermi surface at X or Y in the band structure (inset Fig. 6). Note the discrepancy 
in the empirical i^-value between Figures 5 and 6. A large hidden negative contribution 
to the empirical K can be assumed from the nearest-neighbor Coulomb repulsion V. In a 
superconducting context, Coulomb interactions are reduced to a lower value V* by off-shell 
scattering, while there is no such effect for the static CDW. Hence the CDW effective K 
should be less attractive, with the two K- values differing by roughly V — V* . 

The CDW-induced modification to the energy band structure can be seen by looking at 
the poles in the propagator G(k, n) lying at low energy. Suppose that and ek+Q lie near 
the Fermi level (k and k + Q are 2 points related by nesting), when ek-Q will typically lie 
far away so only the term in ek+Q in $ diverges. Then the poles are approximately at 

iu n = e F ± ^x Q ^k,k+Q -> e F ± xq, (22) 
12 
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FIG. 6: Tcdw (blue circles) calculated from stability boundary of CDW relative to uniform so- 
lution. T c (magenta triangles) and a (green triangles) calculated in presence of CDW from LG 
equation. Inset shows nesting wavevector Q at point X = (vr,0) in BZ. Parameters as Fig. 5, 
except K = 0.23 eV in Tcdw calculation (see text). 



where £k,k+Q ~ ±2. An identical argument can be made for the case of scattering from k to 
k — Q. Hence we see that xq is the gap opened up by the breakdown in long range order 
caused by the CDW - it is a form of pseudogap. The factor £k,k+Q, which is approximately 
cos k x — cos k y at the small Q of interest, shows the presence of a d-wave symmetry factor in 
this gap. Because of nesting factors, the gap is not symmetric in X vs. Y, but full symmetry 
in fc-space will reappear on averaging over spatial domains Q = (Q, 0), or = (0, Q). 

We tentatively identify the CDW gap as the ci-wave pseudogap seen in the HTS materials, 
and identify Tcdw with T*. Experimental studies show a highly inhomogeneous spatial 
variation of the gap Q], apparently controlled by the oxygen dopant distribution, which 
we interpret as a strongly impurity-pinned CDW. The regions where the gap is maximal 
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(75 meV - our xq — 60 meV), but with little sign of superconducting coherence peaks, are 
interpreted as the peaks in the CDW amplitude. The length scale ~ 40 A observed is similar 
to that predicted from Q. 

We need to estimate the interaction between the two order parameters, superconducting 
and CDW. We calculate the CDW order parameter on the basis of the normal state, as- 
suming that always Tcdw > T c (T* > T c ), for which there is supporting evidence The 
superconducting phase - due to its small coherence length £ ~ 2 nm - can coexist with the 



competing CDW phase |46|. We shall here describe the inhomogeneous two-order param- 
eter coexistence phase in an approximate manner using a form of Landau- Ginzburg (LG) 
approach. The approach does not take into account local distortion of the superconducting 
order parameter by the CDW (local distortions average out, leaving perfect ci-wave, over a 
CDW wavelength). 

The LG expression for the superconducting free energy takes the form 

F = \ (VA) 2 + b -(T- T c °) A 2 + 1a 2 / (x 2 (x)) , (23) 
where x ( x ) is the CDW order parameter 

X (x) = X Q cos (Q.x) ; Q ^ 0, X = Xo/V2; Q = 0, (24) 

A(x) is the superconducting gap, a and b and are LG parameters, T c ° is the transition 
temperature in the absence of the CDW order parameter, and / is a coupling function 
between the two order parameters. 

In the absence of CDW order the coherence length is given by the standard formula 



£lg = V a /bT® = 0.739£g CS ; the BCS coherence length ^cs * s assumed to follow BCS 
scaling relative to the known coherence length of the HTS. The function / is determined 
from solving the gap equation as a function of a spatially uniform CDW amplitude xo, the 
input being the splitting t x — t y = Xo i n nearest-neighbor hopping integral. To limit order 
parameter distortion, a cutoff on the maximum x is inserted, which however does not affect 
the conclusions. 

Differentiating the free energy w.r.t. A we get the Schrodinger-like equation (SE) for the 
gap 

-aV 2 A + V(x)A = eA, (25) 



14 



where V(x) = / (x 2 (x)) — 6T C °, e = — bT. The lowest energy e solution represents the 
highest global transition temperature T c = T. The procedure is then to (1) solve for the 
CDW amplitude and wavevector, (2) calculate the effect on T c of a given magnitude of 
uniform order parameter, and then (3) solve the SE for the transition temperature. 

The results are shown in Fig. 6. We see that there is a significant reduction in T c coming 
from the effect of the CDW. The isotope shift still has the same qualitative behavior, but 
with less variation - because xo provides an additional perturbation competing with the 
effect of chemical potential. The STM results^?! might be more compatible with a 2D 
CDW than the ID CDW assumed in this paper, and indeed tentative results for a 2D CDW 
show much less T c -reduction, and enhanced isotope shift variation, but in the absence of 
pinning the ID CDW is found to be more stable than the 2D one. 

In the strong coupling limit, where the effect of K is larger than the effect of t, preliminary 
studies suggest that the FBM can produce a striped phase and a transition to an insulating 
state. Combined with U, a true Mott transition involving magnetism can result, leading to 
an RVB-like state. 

We believe that this work constitutes a breakthrough in finding a natural phonon-based 
mechanism capable of generating a (i-wave superconducting gap, a c?-wave pseudogap, and 
giving the type of behavior of T c , T*, and isotope shift similar to that observed experimen- 
tally. The expected inhomogeneity of the order parameter is also observed Quantitative 
comparisons should bear in mind that the current picture is a highly simplified single band, 
single interaction parameter one. A future program of work comprises a rather long list, 
including first principles substantiation of the anharmonic oxygen potential and coupling to 
the electronic degrees of freedom, inclusion of the oxygen 2p-bands, calculation of spectro- 
scopic properties such as the superconductivity- induced shift in vibrator frequency, and the 
phonon peaks seen in tunneling spectroscopy of the superconducting state, together with 
understanding the fundamental energy balance - possibly unconventional - underlying HTS. 
A better understanding of the CDW order and its interaction with superconductivity should 
be obtained from a more detailed coherent, quantum-mechanical, theoretical approach. 



[1] Leggett, A.J. What do we know about high Tc ? Nature Physics 2, 134-136 (2006). 

[2] Bonn, D. A. Are high-temperature superconductors exotic? Nature Physics 2, 159-168 (2006), 

15 



and the references therein. 
[3] Tsuei, C.C., & Kirtley, J.R. Pairing symmetry in cuprate superconductors. Rev. Mod. Phys. 

72, 969 (2000), and references therein. 
[4] Norman, M. R. et al., Destruction of the Fermi Surface in high-T c superconductors. Nature 

392, 157-160 (1998). 

[5] Moler, K. A., et al., Magnetic field dependence of the density of states of YBa2Cus06.95 as 
determined from the specific heat. Phys. Rev. Lett. 73, 2744-2747 (1994). 

[6] Lee, P.A., Localized states in a d- wave superconductor. Phys. Rev. Lett. 71, 1887-1890 (1993). 

[7] Taillefer, L., Lussier, B., Gagnon, R., Behnia, K. & Aubin, H. Universal heat conduction in 
YBa 2 Cu 3 6 .9 . Phys. Rev. Lett. 79, 483-486 (1997). 

[8] Chiao, M. et al., Low-energy quasiparticles in cuprate superconductors: A quantitative anal- 
ysis. Phys. Rev. B. 62, 3554-3558 (2000). 

[9] Proust, C, et al., Heat transport in a strongly overdoped cuprate: Fermi liquid and a pure 
d-wave BCS superconductor. Phys. Rev. Lett. 89, 147003 1-4 (2002). 
[10] Sutherland, M., et al. Thermal conductivity across the phase diagram of cuprates: low-energy 
quasiparticles and doping dependence of the superconducting gap. Phys. Rev. B. 67, 174520 
1-11 (2003). 

[11] Bel, R., et al. Test of the Wiedemann- Franz law in an optimally doped cuprate. Phys. Rev. 
Lett. 92, 177003 1-4 (2004). 

[12] D. Achkir and M. Poirier,D. A. Bonn, Ruixing Liang, and W. N. Hardy, Temperature depen- 
dence of the in-plane penetration depth of YBa2Cu306.95 and YBa2(Cuo.9985Zno.ooi5)306.95 
crystals from T to T 2 , Phys. Rev. B. 48, 13184-113187 (1993). 

[13] X. J. Zhou, et al., Dichotomy between Nodal and Antinodal Quasiparticles in Underdoped 
La 2 _ x Sr x Cu0 4 . Phys. Rev. Lett. 92, 187001-187005 (2004). 

[14] Anderson, P.W., The resonating valence bond in La2Cu04 and superconductivity. Science 
235, 1196-1198 (1987); Anderson, P.W. et al., The physics behind high-temperature super- 
conducting cuprates: the 'plain vanilla' version of RVB. J. Phys: Condens. Mat 16, R755-R769 
(2004). 

[15] Mourachkine A., The oxygen isotope effect on critical temperature in superconducting copper 

oxides. Supercond. Sci. Techn. 17, 721-723 (2004). 
[16] Pringle, D.J., Williams, G.V.M., and Tallon, J.L., Effect of doping and Impurities on the 

16 



oxygen isotope effect in high-temperature superconducting cuprates. Phys. Rev. B. 62, 12527- 
12533 (2000). 

[17] McElroy, K., et al. Atomic-scale sources and mechanism of nanometer scale electronic disorder 
in Bi2Sr2CaCu208+d. Science 309, 1048-1052 (2005), and references therein. 

[18] Opel, M., et al., Physical origin of the buckling in CUO2; Electron-phonon coupling and Raman 
spectra. Phys. Rev. B 60, 9836-9844 (1999). 

[19] Gallais, Y., Le Tacon, M., Sacuto, A., & Colson, D. Coupling between quasiparticles and a 
bosonic mode in the normal state of HgBa2CuO,4 +( 2 . cond-mat /0506711 v2 (2006). 

[20] Cuk, T., Lu, D.H., Zhou, X.J., Shen, Z.-X., Devereaux, T.P., & Nagaosa, N. A review of 
electron-phonon coupling seen in the high-Tc superconductors by angle-resolved photoemission 
studies (ARPES). Phys. Stat. Sol. (b) 242, 11-29 (2005). 

[21] Cuk, T., et al. Coupling of the Bi 9 phonon to the antinodal electronic states of 
Bi 2 Sr2Cao.92Yo.o8Cu 2 8+(i . Phys. Rev. Lett. 95, 117003 1-4 (2005). 

[22] Gweon, G.-H., et al., An unusual isotope effect in a high-transition-temperature superconduc- 
tor. Nature 430, 187-190 (2004). 

[23] Zhou, X.J.,et al. Multiple bosonic mode coupling in the electron self-energy of 
(Laa^Sr^CuO^ Phys. Rev. Lett. 95, 117001 1-4 (2005). 

[24] Pintschovius, L., et al. Oxygen phonon branches in YBa2Cu30y. Phys. Rev. B 69, 214506 
1-11 (2004). 

[25] Saini, N. L., et al., Evidence for local lattice fluctuations as a response function of the charge 
stripe order in the Lai.48Nd .4Sro.i2Cu0 4 system. Phys. Rev. B 64, 132510 1-4 (2001). 

[26] Lee, J., et al, Vizualizing the influence of nanoscale electron-lattice interactions on high-T c 
superconductivity in Bi2Sr2CaCu20g+«5 (to be published). 

[27] Khasanov, R., Shengelaya, A., Morenzoni, E., Conder, K., Savic, M.I., & Keller, H. The oxygen 
isotope on the in-plane penetration depth in cuprate superconductors. L. Phys.: Condens. 
Matter 16, S4439-S4455 (2004). 

[28] Bussmann- Holder, A., & Bishop, A.R., Anharmonicity-induced multiphonon processes in high- 
temperature superconductors. Phys. Rev. 44, 2853-2856 (1991). 

[29] Bussmann-Holder, A., & Keller, H. Polaron formation as origin of unconventional isotope 
effects in cuprate superconductors. Eur. Phys. J. B 44, 487-490 (2005). and the references 
therein. 



17 



[30] Muller, K.A., From Single to Bi-polarons with Jahn- Teller character and metallic cluster- 
stripes in hole-doped cuprates. To be published in "Treatise on High Temperature Supercon- 
ductivity", edited by J.R. Schrieffer. 

[31] Crespi, Vincent H., & Cohen, Marvin L. Anharmonic phonons and high-temperature super- 
conductivity Phys. Rev. 48, 398-406 (1993). 

[32] Kulic, M.L., Interplay of electron-phonon interaction and strong correlations: the possible 
way to high-temperature superconductivity. Phys. Repts. 338, 1-264 (2000). 

[33] Cappelluti, E. & Pietronero, L., Nonadiabatic superconductivity: The role of van Hove sin- 
gularities. Phys. Rev. B. 53, 932-944 (1996). 

[34] Sakai, T., Poilblanc, D., & Scalapino, D., Hole pairing and phonon dynamics in generalized 
two-dimensional t - J Holstein Models. Phys. Rev. B. 55, 8445-8451 (1997). 

[35] Seidel, A., Lin, H.-H., & Lee, D.-H., Phonons in Hubbard ladders studied within the framework 
of the one-loop renormalization group. Phys. Rev. B. 71, 220501-220505 (2005). 

[36] Schuttler, H.-B., & Pao, C.-H., Isotope effect in d-wave superconductors. Phys. Rev. Lett. 75, 
4504-4507 (1995). 

[37] Lines, M.E., and Glass, A.M., Principles and Applications of Ferroelectrics and related Ma- 
terials, Clarendon Press, Oxford UK (1979), and references therein. 

[38] Biichner, B., Breuer, M., Freimuth, A. & Kampf, A. P., Critical buckling for the disappear- 
ance of superconductivity in rare-earth-doped La2- 2: Sr a: Cu04. Phys. Rev. Lett. 73, 1841-1844 
(1994); Chmaissem, O. et al., Scaling of transition temperature and Cu02 plane buckling in 
a high-temperature superconductor. Nature 397, 45-48 (1999). 

[39] Read, N., and Newns, D.M., On the Solution to the Coqblin-Schrieffer Hamiltonian by the 
Large-N expansion Technique. J. Phys. C. 16, 3273-3295 (1983); Bickers, N.E., Review of 
techniques in the large-N expansion for dilute magnetic alloys. Rev. Mod. Phys. 59, 845-939 
(1987), and references therein. 

[40] Coleman, S., Aspects of Symmetry. Cambridge, New York, (1985). 

[41] Lavagna, M., Millis, A.J. & Lee, P.A., d -wave superconductivity in the large-degeneracy 
limit of the Anderson lattice. Phys. Rev. Lett. 58, 266-269 (1987); Kotliar, G. k Liu, J., 
Superconducting Instabilities in the Large-U Limit of a Generalized Hubbard Model. Phys. 
Rev. Lett. 61, 1784-1787 (1988); Sudb0, A., and Houghton, A., Charge fluctuations, Spin 
fluctuations and superconductivity in a Cu02 sheet. Phys. Rev. B. 42, 4105-4116 (1990). 

18 



[42] Serene, J. W., & Hess, D. W., Quasiparticle properties of the two-dimensional Hubbard model 
in a propagator-renormalized fluctuation-exchange approximation. Phys. Rev. B. 44, 3391- 
3394 (1991). 

[43] Presland, M. R., Tallon, J. L., Buckley, R. G., Liu R. S. , k Flower, N. E., General trends 
in oxygen stoichiometry effects on Tc in Bi and Tl superconductors. Physica (Amsterdam) 
176C, 95-105 (1991). 

[44] Griiner, G., The Dynamics of Charge Density Waves. Rev. Mod. Phys. 60, 1129-1138 (1988), 
and references therein. 

[45] Shibauchi, T., Krusin-Elbaum, L., Li, Ming, Maley, M. P., & Kes P. H., Closing the Pseudogap 
by Zeeman Splitting in Bi2Sr2CaCu2C , 8+ ?/ at High Magnetic Fields. Phys. Rev. Lett. 86, 5763- 
5767 (2001); Timusk, T. & Statt, B., The pseudogap in high-temperature superconductors. 
Rept. Progr. Phys. 62, 61-82 (1999). 

[46] Machide, K., Koyama, T., & Matsubara, T., Theory of charge-density- wave superconductors. 
Phys. Rev. B 23, 99-105 (1981). 

Acknowledgments 

Both authors contributed equally to this work. Correspondence should be addressed to 
D.M.N at dennisn@us.ibm.com 



19 



